Multi-stability in an optomechanical system with two-component Bose-Einstein condensate 
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We investigate a system consisting of a two-component Bose-Einstein condensate interacting dispersively 
with a Fabry-Perot optical cavity where the two components of the condensate are resonantly coupled to each 
other by another classical field. The key feature of this system is that the atomic motional degrees of freedom and 
the internal pseudo-spin degrees of freedom are coupled to the cavity field simultaneously, hence an effective 
spin-orbital coupling within the condensate is induced by the cavity. The interplay among the atomic center- 
of-mass motion, the atomic collective spin and the cavity field leads to a strong nonlinearity, resulting in multi- 
stable behavior in both matter wave and light wave at the few-photon level. 



PACS numbers: 03.75.Mn, 03.75.Kk, 42.50.Pq, 42.65.-k 

In the past few years, there has been a great surge of interest 
in the nonlinear phenomena associated with the so called op- 
tomechanical systems, which are realized by coupling a me- 
chanical oscillator to an electromagnetic field in a cavity 0]- 
7]. The exploration of these systems has led to many excit- 
ing developments, including self-sustained oscillations J^,|^], 
bistability HQxJ, and optomechanical chaos fTHl . Besides 
these linear optomechanical coupling systems, devices with 
nonlinear optomechanical coupling have been studied in some 
very recent experimental works lfl2i [l3ll as well. 

On the other hand, more recent studies on the cavity quan- 
tum electrodynamics (QED) with an ensemble of ultracold 
atoms, both bosonic Il4 - fl8ll and fermionic flffll . give rise to 
a new platform of cavity optomechanics. In this new regime 
of cavity QED, a cavity field at the level of few or even a sin- 
gle photon can significantly affect the collective motion of the 
whole atomic samples. This allows us to study the nonlinear 
dynamics of ultracold atomic gases in a new domain. 

Among all these nonlinear phenomenon in a cavity op- 
tomechanical system, the bistable behavior is one of the fo- 
cuses of research interest. The optical bistability in optome- 
chanical systems has been studied both in theory 112011 and 
in experiment I2H l22tl . Strong matter wave bistability has 
also been investigated in a spinor Bose-Einstein condensate 
(BEQ ll23ll . In this work, we propose a scheme to exploit the 
multistable behavior in a two-component BEC coupled to a 
Fabry-Perot cavity. Here the two components of the conden- 
sate are coupled by another classical optical field, hence real- 
izing a pseudo-spin half system. The cavity supports a single- 
mode standing wave optical field, which interacts with atoms 
dispersively. When the coupling strengths of the two atomic 
spin components and the optical field are different, the cav- 
ity field will then couple to both the external center-of-mass 
and the internal spin degrees of freedom of the condensate. 
Whereas in previous studies of BEC-cavitysystem, the cavity 
field couples either to the center-of-mass jl4t4l^l20l-[22ll or to 
the spin degrees of freedom |[23ll . but not to both l24ll 7As we 
shall demonstrate, the nonlinear coupling among the external 
and internal states of the condensate and the cavity photons 
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FIG. 1: (Color Online) (a) Schematic diagram of the system. The 
cavity is pumped at rate 77 and the cavity photon decays at rate k. (b) 
Atomic level diagram: The two stable atomic hyperfine ground states 
|a) and \b) are coupled by a classical field with coupling strength J. 
In addition, they interact with the cavity dispersively, resulting to an 
energy shift U a and Ut, respectively. 



leads to multi-stability in both light wave and matter wave. 

In our model, as depicted schematically in Fig. 1(a), we 
consider a quasi-one-dimensional BEC of N atoms in two sta- 
ble hyperfine spin ground states \a) and \b) trapped in an high- 
finesse Fabry-Perot cavity. The two spin states are coupled to 
each other by a classical light field (This can be, for example, 
either a two-photon Raman field or an Radio-Frequency field). 
The BEC and the cavity are in the strong coupling regime of 
cavity QED, that is, the maximum coupling strength between 
a single atom and a single intra-cavity photon, g a and gb, are 
larger than both the amplitude decay rate of the intra-cavity 
field k and that of the atomic excited state. The coupled dy- 
namics of the BEC and the cavity field is driven by continu- 
ously applying a weak pump laser field with frequency a> p and 
amplitude r\ along the cavity axis, taken to be the z-axis. 

Let c (c 1 ) denote the annihilation (creation) operators for 
the single-mode cavity field, and a (a') and b (b^) the an- 
nihilation (creation) field operators for the atoms in state \a) 
and \b), respectively. Then the Hamiltonian in a frame rotat- 
ing at the pump frequency oj p can be written as (to focus on 
the atom-photon interaction, we have neglected the atom-atom 
collisions) 
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Here L is the length of the cavity and 6 C = oj p - a> c is the 
cavity-pump detuning. The second line of Eq. (fl} represents 
the coupling between the two spin states by the classical field, 
with / being the coupling strength (without loss of generality, 
we take J to be real and negative), t/,- = g 2 J{u) p — w,) (i = 
a,b) characterizes the cavity field-induced energy shift of 
the atomic spin states, with g, being the resonant coupling 
strength between the atom and the cavity field, and w, the tran- 
sition frequency for the two atomic states, respectively. We 
assume that the pump-atom detuning is large enough so that 
the atomic upper level can be adiabatically eliminated and the 
interaction between the cavity photon and the atom is essen- 
tially of dispersive nature. 

The standing-wave cavity mode couples different momen- 
tum modes of the condensate separated by an integer mul- 
tiple of 2ftk, with k being the wavenumber of the cav- 
ity photon. Recent experiments l2ltl suggest that, when 
cavity photon number is not large, the atomic momen- 
tum modes interacting significantly with the cavity field are 
those with momenta and +2M. Neglecting all higher 
order modes, we can expand the atomic field operators in 
fc-space as a(z) = (ao + yf2a 2 cos 2kz) / V£ and b(z) = 

(b + V2fc 2 cos 2kzj I yfL. Substituting them into Eq. (HJ leads 
to the Hamiltonian under this few-mode approximation: 

— = a>(a 2 a 2 + b 2 b 2 ) + J(a Q bo + <* 2 bi) + J(b ao + b 2 a 2 ) 
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where a> = 4ftk 2 /2m is the photon recoil frequency. 

Now we can draw an analogy with the optomechanical sys- 
tem by defining two harmonic oscillator modes with displace- 
ment^ = {a Q a 2 +a' 2 ao)l ^2N a andXt, = (btb 2 +blbo)/ yJ2N b . 
Their corresponding conjugate variables are P a = -i(a' ci2 - 

a 2 ao)/ ^2N a and P b - —i{btb 2 — b' 2 bo)l \l2Nb~, respectively. 
N s denotes the number of atoms in spin state \s = a,b). Un- 
der the condition that most of the atoms are still in the zero- 
momentum mode, which is an excellent approximation for the 
parameters of current experiments 12 ill , it is straightforward to 
verify that these variables satisfy the commutation relations 
[Xj, Pj] - i(N jo - Nj 2 )/Nj - i (j = a, b), and all other com- 
mutaters vanish. 

In the following we adopt a mean-field treatment by re- 
placing the operators flo(^o) and a 2 (b 2 ) with their correspond- 
ing C numbers ao = y/N ao e~' ea o (j3o = jNh~ e~ ab o) and 
C2 = ^N a ^e~' 6n - (J3 2 = yTV^e - ' 6 ^). Now we can define an- 
other pair of conjugate variables as the atomic collective spin 



M — N a — Nb and the relative phase 9 — 8 ao - 8b . Obviously 
we have [X h M] = [P u M] = and [X t , 6] = [P if 0] = (i = 
a, b). In other words, we have three pairs of independent con- 
jugate variables: (X a , P a ), (X b , P b ) and (M, 6). The first two 
pairs represent the oscillator modes and originate from the 
center-of-mass motion of the condensate, while the last pair 
originate from the internal motion of the condensate. We can 
then rewrite the Hamiltonian under the few-mode approxima- 
tion in terms of these three pairs of conjugate variables as 

-(X 2 + P 2 a +X 2 b + P 2 b ) + JyN 2 - M 2 cos6 



ft 



1 

H — 

2 



(U a - U b )M/2 + U, 



aX a sji 



X b J(N 



(N + AO/2 



c'c + irj (c 1 - c) , (3) 



where A f = 26 c - (U a + U b )N/2. 

To proceed further, we treat the leakage of cavity pho- 
tons phenomenologically by introducing a cavity decay rate 
k whose typical value (~ 1 MHz) is much larger than u 
and 7, under which condition we can assume that the cavity 
field always follows adiabatically the atomic dynamics. From 
ihc = [c, H] = 0, we obtain the mean intracavity photon num- 
ber N c as 

... n 2 
<c t c) = 
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U m m + U a x a 



+ U b x b 
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(4) 

Here we have adopted the normalized variables x a = X a j y[N, 
Xb = Xb/ V^V, m = M/N and the dimensionless parameters 
77 = Tj/K, U a = U a N/K, U h = U b N/K, U,„ = (U a - U b )/2, 
A c = A c /k. For convenience of the following discussion, we 
also define two other dimensionless parameters as oj — cjN/k, 
J = JN/k. 

We can now write down the equations of motion as follow: 
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Pa 


= -U)X a - U a + m)/8 (c 


c), 


(5d) 


Pb 


= -a)Xb - U b V(l -m)/8(c 1 


c), 
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= 2Jm cos 0/Vl -m 2 -\u„ 


+ u a . 


c a /V8(l + m) 



UbXbl V8(l -m)\ {c'c) . 
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where the time derivatives are taken with respect to the dimen- 
sionless time t = Kt. 

Combining Eqs. © and ([5]), one can find the effective 
Hamiltonian H e s as 



-I— = ^r(x l a + P a +x%+p l b ) + Jy\-m 2 cos,6 

TlK 



—rf arctan 
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FIG. 2: (Color Online) Phase diagram for different types of solutions 
with 6 = in the parameter space of ff and A c . Different regions 
are differentiated by their colors. The digits labeled in each region 
denote the numbers of corresponding solutions. The dimensionless 
parameters are taken to be U a = 50, = 200, a) = 400, / = -1. 
The total number of the atoms are set as N = 10 5 . The red dashed 
line corresponds to rj 1 = 7. 
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From Eqs. Q one can see that the atomic dynamics will 
depend on the intra-cavity photon number. Conversely, the 
intracavity photon number is dependent upon the atomic vari- 
ables via Eq. (0). It is this interdependence of the atomic and 
photonic modes that leads to interesting multistability of this 
coupled system, which will be the main result of this work. 

The steady state solutions are obtained by setting the time 
derivatives to zero in Eqs. (01 which yields six coupled nonlin- 
ear algebraic equations. In the case of U a = Ub, i.e., when the 
light shifts induced by the cavity photon are spin-independent, 
it is not difficult to see that m — in the steady state. In 
other words, under this situation, the populations in the two 
spin states are always equal to each other, which effectively 
freezes the atomic spin degrees of freedom. Thus, we will 
always focus on the case where U„ + Ub- 

We may get some useful information from the phase dia- 
gram identifying different types of phase-dependent solutions. 
It follows from Eq. Q that 9 — or n in the steady state. We 
will focus only on the 6 = branch, which for J < repre- 
sents the lower energy branch. In the parameter space of ff 
and A c , the number of steady-state solution are illustrated in 
Fig. |2 We can see that, in certain parameter regimes, the num- 
ber of different solutions of the system can be more than one, 
which indicates that multi-stable behavior may be observed. 
By varying t\ and/or A c , one can traverse different solution re- 
gions of the system. So, this coupled system can be easily 
manipulated by tuning the intensity or frequency of the pump 
laser field. 

As an example, we consider the case where the pump inten- 
sity is fixed at fj 1 = 7. By varying the cavity-pump detuning 
Ac, the equilibrium properties of the system are changed, as 



FIG. 3: (Color Online) Mean intracavity photon number N c (upper 
panel) and normalized collective spin m (lower panel) versus cavity 
pump detuning A c for the steady-state solutions with rj 1 = 7, corre- 
sponding to the red-dashed line in Fig. 2. The branches represented 
by the blue (darker) lines correspond to the dynamically stable so- 
lutions (the solid parts denote the ground state and the dashed parts 
denote stable but not the ground state), and the ones represent by the 
red (lighter) dot-dashed lines correspond to the dynamically unstable 
solutions. 

shown by the red-dashed line in Fig. The corresponding so- 
lutions are derived and the typical results are shown in Fig. [3] 
From these plots, one can see that, both the cavity field and 
the atomic spin population exhibit multistable behavior. For 
certain values of detuning A c , it supports three or five steady- 
state solutions. A standard linear stability analysis shows that 
in the region with three solutions, two of them are dynamically 
stable and in the region with five solutions, three of them are 
dynamically stable. Hence these represent bi- and tri-stable 
regimes, respectively. In these multi-stable regimes, we cal- 
culate the energies of the stable states according to Eq. ©, 
from which we identify the ground-state solution which are 
represented by the solid lines in Fig. [5] As can be seen, the 
ground state jumps from one branch to another at certain crit- 
ical values as A c is scanned. These critical points correspond 
therefore to first-order transitions in this system. 

We remark that optical tristable behavior has been dis- 
cussed in theory several decades ago l26ll and have been ob- 
served in various systems l27lj3"oll . Here we report a new plat- 
form where simultaneous bi- or tri-stability in optical and mat- 
ter waves can be observed. We emphasize that, in the regime 
of weak cavity field such that the few-mode approximation 
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is valid, the existence of tristable regime requires both exter- 
nal center-of-mass and the internal spin degrees of freedom of 
the condensate to be present. By contrast, in previous studies 
where one or the other of these de gree s of freedom is absent, 
only bistable behavior is observed II 1 44 171 1 1 9142311 . 

Note also that the presence of several degrees of freedom 
does not guarantee instability. For example, optomechanical 
systems with multiple de gree s of freedom have been studied 
in several works recently 11251 l3ll l32tl . with no tristability be- 
ing reported. There is a common feature in the systems stud- 
ied in these works, that is each mechanical degrees of freedom 
is only coupled to the cavity mode, and there is no direct cou- 
pling among themselves. In fact, this corresponds to setting 
J — in our model, so that the collective spin M becomes 
a conserved quantity and hence effectively freezes the spin 
degrees of freedom. Then, in principle, through a canonical 
transformation, one can always rewrite the Hamiltonian (see 
Eq. (fj) for example) as a sum of several independent subsys- 
tems, each of which can be regarded as an oscillator coupled 
with a cavity field or just a free oscillator. Obviously, there 
would be no multi-stable behavior other than bistability in 
such a system. 

In contrast, in the model we have considered here, the me- 
chanical degrees of freedom represented by X a and Xi,, the 
spin degrees of freedom represented by M and the cavity field 
are all coupled simultaneously, as can be seen from Eq. (0). It 
is this nonlinear coupling that makes tristability possible. 

In summary, we have studied the interaction of a two- 



component BEC with a standing-wave cavity field, where the 
two components are coupled by another classical optical field. 
We show that this coupled cavity-BEC system can display si- 
multaneously optical multi-stability at the few-photon level 
and matter-wave multi-stability involving a whole condensate 
with a macroscopic number of atoms. This highly control- 
lable optical and collective spin multi-stability can be very 
useful both in exploring fundamental quantum physics such 
as understanding decoherence in a macroscopic system and 
in applications such as building switches and logical gates 
for quantum information processing. Our study also opens 
up possibilities to explore nonlinear dynamical effects such 
as chaos l33ll and bifurcation |33] in an optomechanical sys- 
tem at the regime of few-photon level. In this work, we 
have adopted a few-mode approximation. As a self-consistent 
check, we have verified that for the parameters we used, 
the zero-momentum atomic population always exceeds 90%, 
which should make the approximation valid. In the future, it 
will be instructive to numerically study the validity regime of 
the few-mode approximation and explore the potentially inter- 
esting physics beyond the approximation. Another interesting 
revenue of research is to explore the novel physics induced by 
the coupling between the motional and spin degrees of free- 
dom within the condensate. 
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